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Abstract 

The integrable, (1+1) Gross-Pitaevskii (GP-) equation with hermitian property is extended to chaotic 
behaviour as part of general complex fields within the sl(2,C) algebra for Lax pairs. Furthermore, we prove 
the involution property of conserved quantities in the case of GP-type equations with an arbitrary exter- 
nal potential. We solve the corresponding zero-curvature condition, following from compatibility of mixed 
partial derivatives with Lax pair matrices, in terms of the 'at) '-operators which are derived from a Cartan- 
Weyl basis of the general sl(n, C) algebra or one of its sub-algebras. A gauge invariance of the Lax pairs 
and its zero-curvature relation is proven so that one can reduce the total Cartan-Weyl basis for the spatial 
Lax matrix to the maximal commuting Cartan sub-algebra. This allows to attest the involution property of 
conserved quantities under very general conditions apart from the well-known, classical Y-matrix approach. 
We generalize the approach of Lax pair matrices to arbitrary spacetime dimensions and conclude for the 
type of nonlinear equations from the structure constants of the underlying algebra. One can also calculate 
conserved quantities from loops within the (N-l) dimensional base space and the mapping to the manifold of 
the general SL(n, C) group or a sub-group, provided that the resulting fibre space is of nontrivial homotopic 
kind. This condition is of crucial importance for proving the corresponding involution property and avoids 
possible contraction of the considered loops to trivial point mappings. 
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1 Introduction 

In the study of nonlinear field equations there frequently occur various Lax pairs which allow to construct 
conserved quantities as part of an integrable system. Although it is generally difficult to determine for a 
given nonlinear field equation a corresponding Lax pair, the zero-curvature condition enables deep insight into 
integrability and soliton solutions by the inverse scattering and Backlund transformations for several, physically 
relevant equations, as e.g. the Gross-Pitaevskii (GP)-equation, the Korteweg-de Vries (KdV)-type or the Sine- 
Gordon-like equations pQ-[6]. However, this paper is initiated by a different notion of the integrable (1+1) 
GP-equation which can be summarized by the questions El [9] : 

• Is really any kind of the (1+1) GP-equation (11, ip . e.g. with an external potential V(x,t), integrable in 
the sense of an infinite number of independent, conserved quantities ? (cf. the definition of Liouville 
integrability in Ref. [10] ) 

• Which conclusions can be achieved from the corresponding Lax pair with respect to the involution of the 
conserved quantities and the classical r-matrix ? 

i(d t tp(x,t)) = -(d x d x iP(x,t)) -2V(x,t)i;(x,t) -2\ij(x,t)\ 2 ij(x,t) . (1.1) 

In particular we examine the (1+1) nonlinear GP-equation which has been intensively computed with an 
external potential, a so-called kick- or delta-function potential, for chaotic behaviour on a periodic circle 
ip(x,t) =4>(x + 2Tr,t) [HQ2] 

+oo 

6(t-n) i>(x,t)+g\^(x,t)\ 2 ^(x,t) . (1.2) 

n=— oo 

In advance, we mention that the obvious, hermitian discretization (|1.3j) of the continuous (1+1) GP-equations 
(|l.lll.2p fails to give a corresponding Lax pair whereas the peculiar, discrete, non-hermitian form (jl.4p can be 
obtained from a Lax pair [U EJ [13] 

x (d t ip(x n ,t)) = - (i>(x n+1 ,t) -2iff(x n ,t) +V>(a?n-l)*)) + (1-3) 

- 2 V(x n , t) ip(x n ,t) - 2 \ip(x n , t)\ 2 ip(x n ,t) ; 

i (d t ij)(x n ,t)) = - ^ (V(x n+ i,t) - 2 ip(x n ,t) + ip(x n -i,t)^ + (1.4) 

- (V(x n+1 ,t) + V(x n ,t)) t/j(x n ,t) - |V>(x n ,t)| 2 (t/j(x n+ i,t) +ip(x n -. 1 ,t)) . 

This suggests the preference of the discrete, non-hermitian kind (jl.4p instead of the discrete, hermitian form 
(jl.3p in order to prove an infinite number of conserved quantities from Lax pairs. 

According to a mathematical view, it is necessary to apply the discrete, non-hermitian form (jl.4p of the 
GP-equation in order to preserve exactly the Liouville integrability from the Lax pairs of the continuous kind 
(ll.lll.2p of the GP-equation. Apart from the involution of the conserved quantities which is proven for more 
general cases in later sections, this can explain the numerically observed chaos of the kicked GP-equation 
which may be caused by a false representation of discrete grid points for the numerical integration despite of 
its hermitian form. Therefore, one should use the discrete kind (|1.4p with non-hermitian property in order to 
keep the Liouville integrability from the Lax pairs for the continuous, kicked GP-equation (jl.lll.2p without 
any chaos. 



4 



1 INTRODUCTION 



From a physical point of view, one has to remark that any physical problem has a time-, length-scale or 
energy-, momentum-scale so that the application of delta- function-like time kicks within the (1+1) GP-equation 
has to be regarded as problematic because these time kicks involve any energies (from highly excited Rydberg- 
atoms to cosmic dimensions) |14j . Furthermore, it is not astonishing that different, discrete realizations of 
continuous nonlinear equations can cause different integrable or chaotic behaviour because of varying numerical 
properties; however, in order to maintain the Liouville integrability from Lax pairs of a equation as the kicked, 
continuous GP-equation (|1.2p . one has simply to introduce a finite time-, length-scale or broadened delta- 
function of time for the external kick-potential and has even more to include "many more" (sufficient), discrete 
spacetime steps in the numerical integration so that one can tract the time development for the particular scale 
of the external potential without artifacts coming from below the physical scale. 

In this paper we take the latter point or physical view and furthermore seek for a chaotic behaviour of the 
(1+1) GP-equations. Chaotic behaviour seems to be completely excluded because one can easily conclude for 
non-chaotic properties from the hermiticity and finite norm in any case of a (physical) external potential on 
a spatial, periodically confined circle. Therefore, any kind of GP-equation only seems to allow for integrable 
properties. However, we can also extend to a possible chaotic behaviour , as one considers the su(2)- or sp(2, R)- 
Lax pairs (for an attractive or repulsive interaction) as sub-algebras of the general sl(2, C) algebra and Lax 
pairs. We specify the derivation of a Lax pair for an arbitrary external potential V(x, t) of the general sl(2, C) 
case so that this Lax pair can be used for the construction and involution of the conserved quantities following 
from the monodromy matrix. 

In section [3] we further analyze the point that the construction of a Lax pair, as e.g. for the GP-equation, 
appears to be accidental as one chooses two matrix potentials £(x,t), 9Jl(x,t) for a zero-curvature condition 
which finally determines the nonlinear equations. In appendix [A] we therefore demonstrate the reduction of the 
very general n x n, gl(n,C) Lax pair algebra to its nontrivial, traceless parts sl(n,C); furthermore, we prove 
a gauge invariance of the Lax pairs and their zero-curvature condition in section 13.21 provided that a Maurer- 
Cartan equation holds for the gauge matrices. The gauge invariance of the zero-curvature condition underlines 
that a Lax pair should not be assigned to a single, definite nonlinear equation (as e.g. the characteristic 
GP-case), but to a whole set of equivalent nonlinear equations. We can solve the zero-curvature condition for 
general sl(n, C) Lax pairs (or reduced to one of its sub-algebras as su(n), etc.) in terms of the 'at>'- operator 



for the Cartan-Weyl basis of the prevailing generators of a chosen, closed algebra |15j ; this generally attests 
that the spatial- and time-development of (seemingly accidental) nonlinear equations as part of a whole set 
of equivalent equations is confined to the development within group manifolds, either by the general sl(n, C) 
Lax pair algebra or one of its possible sub-algebras as the su(n)-case. According to this property, it is more 
appropriate to classify the whole set of equivalent nonlinear equations as the GP-like case by the underlying, 
chosen sub-algebras C sl(n, C) for the Lax pairs. 

Since the spatial and time development of nonlinear equations from Lax pairs is limited to that within group 
manifolds, it is not astonishing that one can confirm involution properties of the conserved quantities from the 
monodromy matrix under rather general assumptions by performing the possible gauge transformations to a 
diagonal, spatial Lax matrix £(x,i) — > £^(x,t). This gauge transformation allows to detect conditions for the 
general independence of conserved quantities derived from the Lax pair and its zero-curvature relation apart 
from the introduction of classical r-matrices. Nevertheless, we also investigate the calculation of the classical r- 
matrix in the case of spatially-ultralocal Lax matrices £(x, t) from the tensor product within canonical Poisson 
brackets (cf. section I3~2l and appendix iBl) 



[ki H l + ? a (x,t) E% + Q a (x,t) £° ,...]_.. . 



(1.5) 




(1.6) 
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which can be separated into an eigenvalue part and a commutator of the Lax matrices £i(xi, i) = £(xx,t) <g> 1, 
£2(^21 1) = 1 ® £(x2,t) with the ti2-,t2i- matrices (concerning the general definitions of the tensor product 
space and the tensor product within Poisson brackets, etc. , see Refs. [31 [10], appendix A and chap. 2.5, 
respectively) . 

In section [5] we extend the zero-curvature relation of Lax pairs beyond (1+1) dimensions and determine the 
nonlinear field equations which are also confined to the spacetime development within group manifolds by the 
'oO'- operators of the sl(n,C) algebra or one of its sub-algebras. The calculation of conserved quantities and 
their involution as in the (1+1) dimensional case can be straightforwardly generalized for the corresponding Lax 
pairs within arbitrary spacetime dimensions; however, it has to be taken into account the following additional 
point : 

As we specify closed loops within group manifolds for determining independent, conserved quantities, it has to 
be ascertained that these loops, contained within the mapping of n-1 spatial coordinates to the fields within the 
considered group manifold, possess a nontrivial homotopy in order to prevent possible contractions to trivial 
point mappings. 



2 Lax pairs for GP-type equations from the sl(2, C) algebra 

2.1 Lie algebras of su(2)-, sp(2, R)- and sl(2,C)-Lax pairs for integrable and chaotic be- 
haviour 

In this section we follow Refs. [H |2], but additionally emphasize the algebraic sl(2,C) properties of the Lax 
pairs for the GP-equation. One starts out from two traceless 2x2 matrices £ s i(2,C)(Q) k), 97tsl(2,C)(Q> 7', k) 
()2.4|2.5p with complex- valued entries Q(x,t), l P(x,t) for the physical wavefunctions and a complex- valued 
spectral parameter k (12. 3D which is included to obtain general sl(2, C) matrices. In the context of eqs. (I2.1|2.2p . 
the Lax matrices with the physical fields Q(x,t), 7{x,t) appear to be accidental, but separately determine 
the spatial and time evolution of a two component, complex- valued, (non-physical), auxiliary vector field 
E(x,t) = {E 1 (x,t) , E 2 {x,t)) T 

~ % Q k * ) E(x,t) = £ sl( 2 )C )(a,y;fc) sM) ; (2.1) 

E = E(x, t)eC; 7 = 9(x, t) G C ; Q = Q(x, t) e C ; k £ C ; (2.3) 

£ s i( 2 ,c)(Q,P;fc) = -(**;) [h)+y{x,t) (f+) + Q(x,t) (f_) =^esl(2 5 C) ; (2.4) 

£Wd(2,C)(Q,3 , ;fc) = (2ik 2 + t?Q) (f 3 ) + {-2ky-iy x ) (f + ) + {-2kQ + iQ x ) (f_)^esl(2,q; (2.5) 
; ladder operators (t + ), (f_) of Pauli matrices (ti), [t^), (fsj . 

It is the separate evolution of space and time (|2.1l2.2p for the auxiliary field E(x,t) which enforces the 
pairwise equivalence of the mixed and exchanged spacetime derivatives (I2.6|2,7j) for sufficiently continuous, 
(non- multivalued!) functions. After substitution of (d x E) and (dfE) by (|2.1l2,2p . one transfers the equiv- 
alence of mixed and exchanged derivatives of the auxiliary field E(x,t) to the Lax matrices £ s i( 2j c) (Qj^i^Oj 
3^si(2,c)(2) -Pj k) and attains the so-called zero-curvature relation (|2.8p which restricts the physical fields Q(x,t), 
"?(x, t) to a nonlinear field equation 



zero-curvature condition 
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2 LAX PAIRS FOR GP-TYPE EQUATIONS FROM THE SL(2, C) ALGEBRA 



(d t d x E(x,t)) = (5 t £ sl(2iC) (Q,?;fe)) E(x,t)+S^ 2 ,C)W> k ) ( 5 * 3 (^*)) ( 2 - 6 ) 

= (^Aa(2,C) (Q, JP; *)) t) + £ sl(2jC) (Q, 3>; fc) 9JI sl(2jC) (Q, IP; fc) S(ar, t) ; 

(MS(ar, t)) = ( Wi( 2 ,C) (Q, CP; fc)) S(x, t) + 9K sl(2;C) (Q, 0>; k) (d x E(x, t)) (2.7) 

= {d x m sm (Q, 0>; fc)) S(x, t) + OT sl(2iC ) (Q, k) £ sl(2)C) (Q, 0>; fc) S(x, t) ; 

=> o = ($£a( 2> c)(Q, - (Wi(2,c)(Q,y;A:)) + [^1(2,0(2,^^ , aji sl(2i c)(Q, a 5 ; (2.8) 

As we calculate every matrix element of the (2 x 2) zero-curvature condition (|2.9|) . we (seemingly accidental) 
achieve two coupled, nonlinear equations (|2.10|2,lip of Q(x,t), ( P(x,t) from the off-diagonal matrix entries 
of (|2.8p whereas the diagonal part completely adds to zero without a physical meaning as a wave equation. 
Since the trace of a commutator vanishes for finite dimensional matrices, the zero-curvature condition with the 
commutator between £ s z( 2 j(Q, J"; k) and 9# s i( 2j )(Q, 7; k) always remains within the complex- valued, traceless Lie 
algebra generators sl(2, C). (In appendix|A]we demonstrate how to separate a trivial, diagonal, complex- valued 
unit part from the most general gl(n, C) Lax pair ansatz so that one has only to consider the traceless, sl(n, C) 
parts of Lax pairs and their zero-curvature relation for physical field equations which are indeed nonlinear 
and nontrivial.) It is possible to derive a continuity equation (I2.12|) of the complex fields Q(x, t), 7(x,t) from 
the two field equations (|2.1QI2.1ip which allows to conclude for a time-like conserved quantity with complex 
constant Co as we require spatially periodic boundary conditions (|2.10|2.1ip of the physical fields 

(9t£) (9a- 971) [£,2rt|_ 



o y t \ ( i{7Q) x -2ky x -%y xx \ / %{7Q) X -2ky x -2%y 2 Q , 



Qt J \-2kQ x + iQ xx -i(9Q) x J V ~2kQ x + 2i'?Q 2 -i(?Q) x 

o % + % y xx - 2% y 2 q 



Qt - % Q xx + 2% V 1 = ° 

i (d t y) = (d x d x y) - 2 Q J 32 ; y(x = 0,t) = 9(x = 2vr, t) ; (2.10) 

-i(d t Q) = (d x d x Q) - 2 y Q 2 ; Q(x = 0, t) = Q(x = 2vr, t) ; (2.11) 

=> 9 ( (?Q) = -zd x (Q(d x 9)-y(d x Q)) ; (2.12) 

/"27T /"27T 

St / dz Q(a;,t)) = -» / da; Qb(Q(x,£) (fla-IPfot) ) - 7{x,t) (d x Q(x,t) )) = ; (2.13) 

Jo Jo 

r 2ir 

dx (y(x, t) Q(x, t) ) = const. = c GC. (2.14) 

o 

The conserved, complex-valued quantity (|2.14p refers to a generalized norm of Q(x,t), y(x,t), compared to the 
standard GP-equations, but lacks of its positive definiteness. It is the non-positive property of the generalized 
norm (12.14p which will allow us to accomplish a possible chaotic behaviour within sl(2, C) Lax pairs of physical 
fields. The restriction of the sl(2, C) generators £ s i(2,c)> ^sl(2,C) to the sub-algebras su(2) and sp(2,R) leads 
to the well-known, (1+1) GP-equations (|2.15|2.16p for the attractive and repulsive interaction case. As we 
set in (|2.15p (Q(x,t) = —7*(x,t) = ip(x,t)) to the physical wavefunction with periodic boundary conditions 
and with a real spectral parameter k, we acquire anti-hermitian Lax pairs £ SU (2)(Q; ^> 9Jt su ( 2 )(Q, k) whose 
zero-curvature condition determines the attractive case of interaction of the (1+1) GP-equation 

I) : attractive case, su(2) sub-algebra of sl(2, C) (2-15) 
Q(x, t) = -T(x, t) = V>(x, t) G C ; keR; 
%[dtif>) = -[dM)-2\il>\ 2 i/>; i/>(x = 0,t) =V(x = 2vr,i) ; 



2.1 Lie algebras of su(2)-, sp(2, R)- and sl(2,C)-Lax pairs for integrable and chaotic behaviour 
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^)(^k) = zA; J; ^ su(2) (Q,y;fc) = ^ _ 2 ^ + ^ ^^ij; 

-2 S u(2)(Qi 3*; = -(Asu^CQ,?;*;)) 1 ; ^su(2)(Q,y;fc) = -(stn 8tt(2 )(Q > o j ;fc)) t . 

Similarly, one can choose the combination (Q(x,i) = = vp(x,t) ) for the physical wavefunction with a 

real spectral parameter k in order to achieve Lax matrices within the symplectic algebra sp(2,R), having the 
real parameters ^R(ip(x, t) ), Q(ip(x, t) ), k £ R; the zero-curvature relation then results in the GP-equation with 
a repulsive interaction 

II) : repulsive case, sp(2, R) sub-algebra of sl(2, C) (2-16) 
Q(x,t) = +y*(x,t) = ip(x,t) G C ; {% k — G R ; 

= -(aA^) + 2|^| 2 ^ ; ^(z = 0,t) =V(x = 27r,t) ; 

Z sp (2,R)( Q > 7 > k ) - { ^ k ) ' sp(2,M) (Q) j "<) ^ 2 *A^ + i^ 2z/c 2 - * |-0| 2 J' 

o = £ sp ( 2 ,R)(Q,y;A:) (J "o 1 ) + ( ? 'o 1 ) ( £ sp(2,k)(Q^^)) T ; 
o = 9^ S p(2,R) (Q) k) (J V ) + ( ? "o 1 ) (^ p (2,k)(Q,P;A0) t . 

Both cases, the su(2)- and sp(2, R)-sub-algebras for the attractive and repulsive interaction, have the identical 
continuity equation (|2.17p . due to the hermitian property of the corresponding Hamilton operator. The physical 
norm of the periodic fields on a circle is always positive definite and is set to one (or e.g. to a constant particle 
number Nq) for an interpretation of a true probability density. Apart from the normalization of wavefunctions 
at arbitrary time, we can apply the positive norm of constant value to estimate the maximal achievable value for 
the norm of the difference of two wavefunctions at different times. Since the absolute value of any wavefunction 
is restricted by the maximal probability density one at arbitrary times, the maximal achievable difference of 
the absolute values of two wavefunctions at two different time points is bounded by the value 8ir, due to 
the restricted spatial coordinate x G [0, 2tt). Therefore, the hermitian property of Hamilton operators in the 
attractive su(2)- and repulsive sp(2, R)-cases only allows to assign a non-chaotic behaviour to the GP-equations 
(j2. 1512. 16|) . due to the absence of an exponential- (or other type-)divergence at two arbitrary times 

=► d t {r^) = id x (r (d x ?p)-ip(d x r)) ■, (2.17) 

=^ 3 t I dx(i>*(x,t)^(x,t)) = % [ dxd x (i>*(x,t)(d x i>{x,t))-i>{x,t)(d x ^*{x,t))) =0; (2.18) 



2tt 

dx (ip*(x,t) ip(x,t)) = const. = rg £ R ; r := 1 ; (2.19) 



2tt 



dx (ip*(x,t) - ip*(x,t = 0)) (ip(x,t) -if>(x,t = 0)) = (2.20) 
dx (\*p(x,t)\ 2 + \i/)(x,t = 0)| 2 - 2 cos (Z(il>*(x,t)\il>(x,t = 0))) \ip(x,t)\ \t/j(x,t = 0)|) < 

f'llT flit 

< / dx (1 + 1 - 2 • (-1) • 1 • l) = / dx4 = 8ir. 
Jo Jo 

In section 0] we also prove the time-like conservation of traces from arbitrary powers of related monodromy 
matrices (12.21112.23]) which even fulfill the involution property of the Poisson brackets. According to our sep- 
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2 LAX PAIRS FOR GP-TYPE EQUATIONS FROM THE SL(2, C) ALGEBRA 



arate consideration of the three algebraic cases, one can derive conserved quantities from the T s i(2,c)( a; ) *j k) 
matrices for the rather general sl(2,C) case with independent fields Q(x,t), "P(x,t) or for the two special 
su(2)-, sp(2, R)-sub-algebras with the exponential matrices T su ( 2 ) (x, t; k), 3w 2 ,]R) (^j t) k) °f the- spatial Lax 
generators £ su ( 2 )(£, *j k), £ S p(2,K)(£> The inequality (|2.20p . which restricts the wavefunctions to non- 
chaotic behaviour, complies with the infinite number of independently conserved quantities derived from the 
T su ( 2 )(x, t; /c), X S p(2,]R)(x, t; k) matrices for the attractive and repulsive case, respectively (cf. section 0] and 
appendix [B]). According to the inequality (|2.20p and the independently conserved quantities from the traces 
of powers of T su ( 2 )(x, k), Z S p(2,R){ x i t; k), it seems to be impossible to obtain any chaotic behaviour from the 
su(2)-, sp(2,R)- Lax pairs and their zero-curvature conditions. Chaotic behaviour from Lax pairs appears to be 
impossible in particular because the inequality and the independently conserved quantities from T su (2)(x,t; k), 
T sp (2 R) (x, t; k) can be simply generalized to the case with an arbitrary external potential in the GP-equations 
with attractive or repulsive interaction 



%\(2,c)(x, t-,k) = e*p{J Q d £ £ s i(2,c) 



C sl(2,C)^ k) = Tr ( T sl(2,C)(^ = 2tt, t; k)) 1 



(2.21) 



^sl^.C)^'*^) ^sl(2,C){x, t;k) ; 1 s ^ 2 C) (x,t;k) =exp^- d£ £ s i( 2 ,c)(£, *! k )] ! 



%u(2)(x,t; k) = expj / d£ £ su ( 2 ) 



ci:\ 2) (t;k) = Tr 



su(2)( 



( x ' fc ) T su(2) (x, t; k) ; ^ u(2) (x, t; k) = exp{ 



% P {2,M,){x,t]k) 



-1 

1 



~su(2) 

exp 

sp(2,R) 



= 2tt, *;*;))"] ; 

^su(2) 



(2.22) 



sp(2, 



(x,t;k) 



(£,*;*)} C^} m (t;k) = Tr[(% p{m (x = 27r,t;k)) n ] ; (2.23) 
< ^sp(2,R)( a; ) *j k) ; 



-1 

1 



sp(2,l 



) (^) t\ k) 



CXPi ^ ^£ sp{2j ; 



Nevertheless, it is possible to compose the fields Q(x,t), T(x,t) of the rather general sl(2,C) Lax pair algebra 
by the combination (|2.25|2.26p of the physical wavefunctions ip sn (2)( x i V'sp(2,R)( x ) *) so that probability of 
su(2)- and sp(2, R)- sub-algebra matrices with their corresponding wavefunctions as parameters is allowed to 
flow and change within the total sl(2, C) Lax pair algebra; in consequence chaotic behaviour becomes possible 
despite of the formulation in terms of Lax pairs and their zero-curvature relations 



c (G C) 



2t 



dx (T(x,t) Q(x,t)) 



Q(x,t) = Vsu(2)(<M) +A P {2,K)(x,t) ; 
T{x,t) = -^* u{2) (x,t) +ip* sp(m {x,t) 



(2.24) 

(2.25) 
(2.26) 



We illustrate this possibility from the generalized conserved norm (|2.24j) of fields Q(x,t), ( 3 3 (x,t) with the 
complex constant cq. One has to compose the fields Q(x,t), 7{x,t) of sl(2,C) by the wavefunctions ip su (2)( x j't)^ 
ip Bp (2m(x, t) of the attractive and repulsive interaction case of GP-equations according to eqs. (|2.25l2.26p so 
that one can resolve from the real "5R(co)" and imaginary part "3f(co)" of the generalized norm (|2.24p two 
independent conserved relations (|2.29l2.30p in terms of the independent physical wavefunctions ip su (2)( x i 0> 
'4'sp(2,R)( x ^ t). In correspondence to the non-compact property of SL(2,C), we attain unbounded norms of 
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i ; su(2)( x i t), tp S p(2,R)( x i t), which only differ by the constant K(co), and also acquire a restriction for a generalized 
angle between the wavefunction vectors of i/j SVl (2)( x i *)> VVp(2,R)( x > from the imaginary part (as one transfers 
the spatial argument "x" of wavefunctions to a discrete index). The decomposition of Q(x,t), T(x,t) to the 
parameters V'su(2)( a; j ' l Psp(2,R)( x ^) °f the sub-algebras su(2), sp(2,M) within sl(2,C) therefore removes the 
restriction to completely integrable behaviour (even for the cases with arbitrary external potential V(x,t)) 

Co(GC) = I dx (- V* u (2)0,*) +V , s *p(2,R)( :C >*)) (i>Bu(2)( x ,t) + Ap(2,R)(x,t) 
J 

2tt r 2ir 



dx |V>sp(2,R)(M)| 



dx \Au(2)( x ,t)\ +K(c ) + 



— I 



3(co) 



2m 



cfa; 9(^* u(2) (a;,i) ^ sp ( 2 ,R) (a;, i) 



(2.27) 
(2.28) 



2 77 



277 



dx |^sp(2,R)(M)| 
dx %(lp* u{2) (x,t) Vsp(2,R)0M) 



277 



i 2 

dx \ An(2)(x,t)\ + 3?(c ) ; 



(2.29) 
(2.30) 



After replacing the fields Q(x,i), "P(x,t) by (j2.25j2.26j) in the sl(2,C) zero-curvature relation and field equa- 
tions (|2.3ip . we attain two coupled GP-type equations with the attractive case of the su(2)-field ip su (2)( x , t) 
and with the repulsive case of the sp(2, R)-field ip S p(2,R)( x ,t). Aside from the cubic, attractive interaction 
'—2 |V'su(2)( x ) *)| V'su(2)( x ) 0' an d the cubic, repulsive interaction '+2 \ip sp (2,R)( x i t)\ 2 ' ! /'sp(2,R)( :c ) t)i there appear 
the parts '+4 |V ; S p(2,R)( a:; > 01 VV(2)( X ) an d ' — 4 |V'su(2)( a; ) *)| 2 ^sp(2,R) (^j acting as a repulsive and attractive 
external potential, respectively. Moreover, one has the two incoherent terms ^—2ijj* u( ^{x,t) ^ 2 p ( 2 r)( x ' ana - 
'+2 ^ S p*(2,R)( a; ) ^su(2)( x '*)' wn i cn can give rise to a change of norm and integrated, total probability 



-z (<9 t r) 



9 
Q 



2 r T 2 
23>r 2 



« (9tAu(2) ) 
* (dtV>sp(2,R) ) 



{d x dxAu(2)) - 
-(d x d x Ap(2,R)} 



- 2 
+ 2 



^su(2) | 

^ sp (2,: 



2[Vsp(2,M) | J ^su(2) - 2 T/4(2) (<A 



sp(2,l 



>) 2 ; 



2|^su(2)| ) V> S p(2,R) + 2 V s *p(2,R) (^su(2))' 



This becomes obvious as we derive continuity equations with the corresponding source parts 



{dt\Au{2)(x,t)\ 2 ) 
(3 t |V>sp(2,R)(M)| 2 ) 



-2 a.3 



(^(2)0^*) (dxAn(2)(x,t))) -4 |Vfeu(2)0M)| 2 $ (^ 2 p ( 2 ,R) 0, *) ) ! 
(V>sp(2,R)O>0 (^^sp(2,R) (».*)) ) + 4 IV>sp(2,R)(M)| 2 ^(^(2) ( X > *) . 



(2.31) 

(2.32) 
(2.33) 

(2.34) 
(2.35) 



According to the coupled GP-type equations (I2.32)2.33p . one can distinguish three different, physical cases. 
If we set the two wavefunctions ip BU (2)(x,t) > V'sp(2,R)( x > to an identical value ipo(x,t), one just considers 
the trivial free propagation with the kinetic term , —(d x d x il)o(x,t)y (|2.36p . However, we can also take the 
attractive or repulsive case with the wavefunctions ip SVL (2)( x >t) > V , sp(2,R)( x > as the dominant term ipo(x,t) in 
the two, coupled, nonlinear equations (j2.32|2.33j) and can then regard the other ip sp (2,R)(x, t) or V'su(2)( a; j t) wave 
component as the corresponding perturbation Aip(x,t) in order to compute for a possible chaotic behaviour of 
the i/Jsu(2)( x j t) and ip S p(2,R)( x i *) wavefunctions, respectively 

1) h : ^su(2) = ^s P (2,R) = ; A?/; = free propagation ; i (d t tp ) = -(d x O x ip ) 5 (2-36) 
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2) if : VW2) = tpo ; V'sp(2,E) = (perturbation) chaotic behaviour of V"su(2) ! (2.37) 

3) if : V , S p(2,R) = V'o ! V'su(2) = (perturbation) =^> chaotic behaviour of V ; S p(2,R) • (2.38) 

The latter two cases (|2.37I2.38|) as sub-algebras of sl(2,C) have an infinite number of independent, exactly 
conserved quantities which are derived from traces of powers of Tsi(2,c)( a: ^i k) monodromy matrices. These 
exactly conserved quantities from X s i(2,c)( x ' t; k) contain both sub-groups, the SU(2)-group with !T SU ( 2 ) ( x , t; k) 
matrices and the Sp(2, R)-group with the 3w2,r) (x, t; k) matrices, so that one can consider perturbations of 
the dominantly conserved parts, which are given as traces of powers of 1 Z su (2)( x it] k) or X sp (2,R) (x,t;k) with 
corresponding perturbative terms of £ Bp (2,m.)( x > t] k) or £ fm (2)( :c >*; k), respectively. 

In this section we have limited discussion to sl(2,C) Lax pairs and zero-curvature relations with the su(2), 
sp(2,R) sub-algebras of the corresponding integrable, attractive or repulsive GP-equations. However, the 
coupled equations (|2.32|2.33l) of tp su (2)(, x , t) and ip S p(2,R)( x jt) wavefunctions within the non-compact sl(2, C) 
algebra can be generalized to arbitrary algebras sl(n, C) with a chosen sub-algebra su(n) so that one can 
achieve a possible chaotic behaviour for field equations following from sl(n, C) Lax pairs and zero-curvature 
conditions. The rather general sl(n, C) Lax pair has to be decomposed into sub-algebras as the compact 
su(n) case of anti-hermitian Lax pairs and remaining algebra parts so that the norm of wavefunctions is not 
bounded by real, positive constants. This extension to chaotic behaviour, which follows from taking compact 
sub-algebras of the rather general sl(n, C) algebra, is also possible for Lax pairs beyond (1+1) dimensions. 



2.2 Lax pairs of sl(2,C) for an external potential 

The derivations, concerning the integrable and chaotic behaviour in the previous section 12.11 can even be 
transferred to the case with an arbitrary external potential V(x,t) in the attractive su(2) and repulsive 
sp(2,R) interaction case of GP-equations. In the following we construct sl(2,C) Lax matrices £ s i(2,c)(Q> ^5 k), 
SPT s l(2,c) (A, B, C) so that the nonlinear field equations for Q(x, t), 7{x, t) of the zero-curvature relation allow un- 
der restriction to ^ S u(2)( x > 0> ' ! /'sp(2,R)( x ) t) fields and sub-algebras the GP-equations with the external potential 
V(x,t). It suffices to derive the nonlinear equations 

»JP t = 7 XX -2(Q'P-V(x,t))y ; (2.39) 
-tQt = Q xx -2(9Q-V(x,t))Q, (2.40) 

so that the whole discussion of section [2~TT1 can be conveyed to the case with an external potential, especially the 
conclusions for integrable and chaotic behaviour. We point out that the spatial Lax matrix £ s i(2,c)(Q> k) does 
not change aside from a possible spacetime dependence of the complex spectral parameter k — ► k(x,t) G C. 
Therefore, the conserved quantities, following from traces of powers of the monodromy matrix % B U2 C){ X > ^ k), 
are not essentially altered and the canonical variables Q(x,t), y{x,t) within the Poisson brackets do not cause 
a different r-matrix which specifies the involution of the independently conserved quantities. (In the sequel one 
is even allowed to choose a constant spectral parameter k G C, concerning the involution property, because the 
complex-valued, auxiliary field W can always be adapted in such a manner so that an unwanted appearance 
of V x (x,t) can be absorbed by suitable dependence of W x in (|2.59|) .) 

In analogy to eqs. (12. 1112. 5j) . we start out from traceless matrices £ s i(2,c)(2> ^) k), iJJl s u2,c)(A,B,C) where 
the latter matrix for the time evolution of the auxiliary field H(x, t) is determined by an ansatz with the three 
complex-valued fields A^, Q; k), B(7, Q; k), C(9, Q; k) 

\ E(x, t) = £ sl(2)C) (Q, 9; k) S(x, t) (2.41) 




2.2 Lax pairs of sl(2, C) for an external potential 
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s ' = (c&%k) Tw)) s(M)= ^ IAS ' c)3(, '' ); (2 ' 42) 

S = E(x, t) G C ; 5 = ^(x, i)eC; Q = Q(x, i)eC; fe = k(x, t) G C ; (2.43) 

A = Q; fc) G C ; B = B{y, Q; fc) G C ; C = C(0>, Q; fc) G C ; (2.44) 

£si(2,c)(Q,y;fc) = (fs)+0>(x,t) (f+) + Q0r,t) (f_) =«Esl(2,C); (2.45) 

^si(2,c) (A, B,C) = A(T, Q; fc) (f 3 ) + 5(y, Q; fc) (f+) + C(3>, Q; fc) (f_) =^G sl(2, C) . (2.46) 

The zero-curvature condition (|2.47p results into the field equations (|2.49l2.50p and additionally into the spec- 
tral parameter equation (|2.51|) which follow from the off-diagonal and diagonal matrix entries in f|2.48[) . A 
modified ansatz (|2.52M2.54"j) for the coefficients A(7, Q; k), B{y, Q; k), C(T, Q; k), similar to the case in previous 
section |2~T| leads to the equations H2.57j2.58p of the fields Q(x,t), < ?(x,t). Aside from the external potential in 
the coefficient A(1P, Q; k) (|2.52p . the ansatz (|2. 52H2.54"j) is supplied with a sufficient number of complex- valued 
fields W(5 > , Q; k; V), Y^F, Q; k; V), Z(y, Q; k; V) (|2.55|) which have to be determined in such a manner that the 
nonlinear field equations (|2.57l2.58p of Q(x, t), { J'(x, t) only change by the external potential V(x, t). These con- 
ditions (cf. the two braces in eqs. (|2.57|2.58j) ) specify first order partial differential equations of Y x {7, Q; k; V), 
Z x {7, Q; k; V) whereas the additional equation (|2.59p of a (possibly chosen spacetime dependent !) spectral 
parameter k(x, t) results into a differential equation with the partial derivative W a; (3 5 , Q; k; V) 



zero-curvature condition 

= (d t £ sl(2)C) (Q, ?;*;)) - (d x M Bl{2jC) (A, B , C)) + [£ s i (2 , c) (Q, ?; k) , M sl{2tC) (A, B,C) 
-ik t 9 t \ ( A x B x \ ( "PC-QB -27 A-2ik B 



Q t ik t J \C X -A x J + V 2QA + 2ikC -(9 C - Q B) J= 0; 

y t = B x + 2 9 A + 2i k B ■ 

Q t = C x -2QA-2ikC; 

-ik t = 4-TC+QB; 

A = 17Q + 21 k 2 -iV(x,t) + W ; 

B = -i r S> x -2k r S> + Y ; 

C = iQ x -2kQ + Z ; 

W = W(J>, Q; k; V) G C ; Y = Y(V, Q; k; V) G C ; Z = Z(T, Q; k; V) G C 

k = k(x,t)eC; V = V(x,t)eC; 



z y t = <J> XX - 2 (Q 3> - V) 0> +i Y x - 2i (k x - W) 9 - 2 k Y ; 
-i Q t = Q xx — 2 {y Q — V) Q —i Z x + 2i (k x + W) Q — 2k Z ; 



2.47) 

2.48) 

2.49) 
2.50) 
2.51) 
2.52) 
2.53) 
2.54) 
2.55) 
2.56) 

2.57) 
2.58) 



-ik t = i(2{k 2 ) x -V x ) + W x + YQ-ZT . (2.59) 

We combine the restrictions of the auxiliary fields W{7, Q; k; V), Y{7, Q; k; V), Z(7, Q; k; V) (|235l of the 
modified ansatz with coefficients A{7,Q;k), B{7,Q;k), C(IP, Q;k) from relations (|2.47H2.5"9"|) and require the 
nonlinear field equations (|2.6QI2.6ip with the external potential V(x,t) for the physical fields Q(x,t), '?(x,t). 
The nonlinear field equations (|2.6U|2.61|) only follow from the zero-curvature condition with modified Lax matrix 
9Hb1(2,C)(4 B ' C). provided that the additional fields W(9, Q; k; V), Y(V, Q; k\ V), Z(3>, Q; k; V) (12351) fulfill the 
first order, spatial evolution equation (|2.62p . However, this evolution equation (|2.62p in the spatial variable x 
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can be solved under rather general conditions with the appropriate ordering of exponential step operators and 
initial value iv(x = 0,t) (|2.63|2.64|) . Apart from the homogenous solution with latter field vector ro(x = 0, t), 
one has a particular solution with vector field d(J > , Q; k; V) which also contains the spatial derivative of the 
external potential V(x,t) 



iVt = y xx -2(Q?-V(x,t))?; (2.60) 

-iQ t = Q xx -2(7Q-V(x,t))Q; (2.61) 

(W\ ( -Q 7\(W\ ( i{V x -k t -2(k 2 ) x ) \ 

d x [ Y \ = -27 -2ik \ \ y + 2k x T ; (2.62) 

\ Z ) \2Q 2ik ) \ Z ) \ 2k x Q / 

d x m = W(9,Q;k)t6 + S(7,Q;k;V) ; (2.63) 



W = {W,Y,Z) T ; »(7,Q;k;V) = (i{V x - k t - 2(k 2 ) x ) , 2k x 7, 2k x Q) ; 

^W(^i)} tv(x = 0,t) + J dyexp^J d£ W(^,i)} . (2.64) 

Therefore, one can construct Lax matrices for generalized GP-type equations (|2.60|2.6ip with the additional 
fields W(!P, Q;&; V), Y(T,Q;k;V), ZiT^Q^k^V) (I2.55l2.62ll2.6"4l) in such a manner, that the zero-curvature 
relation (|2.47|2.48[) confines to the cases with an external potential V(x,t). As we have already mentioned 
at the beginning of this section 12.21 one can directly repeat considerations of section 12.11 in order to conclude 
for integrable or a possible chaotic behaviour by regarding the corresponding cases of sl(2, C) Lax pairs with 
sub-algebras su(2) and sp(2,R). 

3 The general n X n Lax pair matrices as sl(n, C) algebras 
3.1 Solution of the zero-curvature condition 

The most general matrices for the Lax pairs £(x,t), 9Jt(x,t) are given by the complex-valued gl(n,C) algebra 
with 2n 2 real parameters; however, due to the occurrence of the commutator [£(x,t) ,9Jt(x,t)]_, one can split 
a trivial, diagonal unity part from this most general ansatz with two complex-valued fields consisting of the 
sum of the diagonal entries from the £(x,t), $Jl(x,t) matrices within gl(n,C) (cf. appendix|A|). Therefore, the 
sl(n,C) algebra with 2 n 2 — 2 real parameters is only taken as the most general ansatz for Lax pairs £,(x,t), 
9Jt(x, t) so that, indeed, field equations remain with nontrivial, nonlinear properties. 

In the following, we assume that the spatial Lax matrix £(x, t) is chosen within a sl(n, C) algebra or 
one of its sub-algebras as e. g. su(n), etc. . Therefore, one can view the zero-curvature relation as a spatial 
evolution equation (13.11) of (d x 9Jl) with the aD-operator [£(£,£) ,...]_ of a closed algebra acting onto Wl(x,t) 
and an inhomogeneity (dt£(x,t)). The spatial evolution equation (|3.ip has a homogenous matrix solution 
Tl hom (x,t) (Hd with 'ini(tial)' matrix Wl in 

■,i(tiai)( x — 0, t) a ^ the spatial origin which develops with spatial 

steps Ax - [£(£,£) , . . .]_ in exponential operators with appropriate spatial ordering from right to left 'exp{. . .}'. 

> 

The action of the aO-operator [£(£, t) , ...]_ within the spatially ordered exponentials has the effect of left 
and right propagation for the initial matrix 9Jtj ra j(x = 0, t) with exponents of ±Ax • £(£,£); this is reminiscent 
of the Heisenberg equation of motion in quantum mechanics or of the total development operator for the 
von-Neumann equation in statistical mechanics 

(d x m) = [£ , sen] _ + (dt£) = [£, + ($£) ; (3.1) 



3.1 Solution of the zero-curvature condition 
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homogenous solution : 9Jth om (x, t) ; 

m*m(x,t) = ^p{j Q ...]-}Tl ini{tial) (x = 0,t) ; (3.2) 

9^hom(3 ; ) i) = exp I Ax £(x, i)} exp { Ax £(x — Ax, t)} exp {Ax£(x - 2Ax,t)} ... x (3.3) 
x... exp {Ax£(2Ax, t)} exp {Ax£(Ax,i)| TX ini (x = 0,t) exp { — Ax £(Ax, t)\ exp { — Ax £(2Ax, t)\ ... x 
x . . . exp { — Ax £(x — 2Ax, t) } exp { — Ax £(x — Ax, t) } exp { — Ax £(x, i) } . 

The particular solution of (|3.ip similarly results as in the case of an ordinary, first order differential equation 
by a variational matrix-ansatz 9T(x,i) (|3.4p instead of the initial matrix 9JTj n j(x = 0, t). Straightforward 
transformations (|3.5I3.6[) lead to the general solution (|3.7p which consists of the homogenous part with matrix 
%Rini(x = 0,t) and the particular part with the inhomogeneity (dt£(y,t) ) 



Variational ansatz for initial matrix 9Jtj n j(x = 0, t) with yi(x,t) 
Wl{x,t) = e^{J^ d£ [£(£, t) , . . . ]-}5Wim0c = 0, t) + ex^{ ^ d£ [£(£,*), ■ ■ .]-}0T(x,t) 



(3.4) 



inserted into zero-curvature condition 



e3q5{ jT d£ [£(£,*), •••]-} (d x m(x, t) ) = (d t £(x, t) ) ; (3.5) 

pgr py 

Sl(x,t) = j dyel${- j dCm,t), ...]_} (d t £(y,i)) ; 



(3.6) 



OJtOz, t) = ex^j £ d£ [£(£,*), ...]-}SKim(a: = 0, t) + ^ e3Ep{ jT ^ [£(£,i), ...]-} (S t £(y, t)) • (3.7) 



The general solution (13. 7j) of the zero-curvature relation demonstrates that it suffices to choose a spatial Lax 
matrix £(£,£) from a closed algebra as the general sl(n,C) case or one of its sub-algebras as su(n) in order 
to determine the matrix 5DT(x, t) or the physical, nonlinear equations wave equations following from it. After 
choosing the spatial matrix £(£, t) and the initial matrix 5D?j n j(x = 0, t) within a closed Lie algebra of a Cartan- 
Weyl basis (|3.8p . one is constrained to the spectral development within the corresponding, closed Lie group, 



due to the action of the aU-operator [£(£, t) , ...]_ of the closed Lie algebra. This points to a classification 
of nonlinear equations by the underlying closed Lie algebra of Lax pairs (as the general sl(n, C) case with its 
sub-algebras as e.g. su(n)) instead by the precise form of the physical, nonlinear equations (as e.g. the typical 
GP-equations). This becomes even more obvious as we can prove a general gauge invariance from a closed 
algebra of Lax pairs so that one has to regard a whole set of Lax pairs for a certain type of physical, nonlinear 
equations. Let us consider the Cartan-Weyl basis (13. 8p for a closed Lie algebra |15| 



[H\W]_ = 0; = N a f>E«+0; (a + ^ 0) ; 

[& ,E a ]_ = a? E a ; [E a , E~ a ] = ay W , ^ ' ' 

so that the spatial Lax matrix has the general, ultralocal form (13.9[) with the fields 4> a (x, t) = (Q a (x, t), T a (x, t) ) 
for the ladder operators E a = (£*" , E9) and with the general, spectral parameters ki(x,t) for the (maximal 
commuting, traceless) Cartan sub-algebra. According to the gauge invariance of Lax pairs (cf. next section 
3.2p . one can transform the general ansatz of £(x,t;k) (|3.9p with the ladder operators E± to the (maximal 

commuting, traceless) Cartan sub-algebra generators H l where the corresponding fields \ { " \x,t;k) (|3.10|) 
follow from the fields <j) a (x,t) and the spectral parameters 

£(x,t;k) = H* h(x,t) + E a </> a (x,t) ; (3.9) 
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E a </> a (x,t) 

£(x,t;k) 
Wl ini ( y x = 0,t) 
E a B ( ^{x,t) 

m(x,t) 

E a b a (x,t) 

m(x,t) 



Tt(x,t) 



E!Q a {x,t) + E%? a (x,t) ; <j> a {x,t) = (Q a {x,t) , 9 a (x,t)) ; 

£&(x,t;k) = H l X[ H } (6 a ; kj) = H l X\ H \x,t; k) ; (by a gauge transformation); (3. 10) 
H i Af ni) (x = 0, t) + £ Q s£ ni )(x = 0, t) ; 

fl* ai{(f)p{x, t); kj) + £ Q b a (<f>p(x, t); kj) = H l ai {x, t) + E a b a (x, t) 
E a _ b-, a (x,t) + E% b+, a (x,t) ; 

J*dy tiw{J*d£>f\S,t) [H\ ...]!} x 

IP {d t \f\y,t))+6(y) (iT Af%,i) + £ Q 
iT ai(^(x, t); kj) + £ Q 6 Q {<f> p (x, t); kj) = H* a^x, t) + E a b a (x, t) 



(3.11) 
(3.12) 

(3.13) 



(3.14) 



W (Af^(x = 0,t)+ I dy(d t Xp(y,t))) + 



(ini) 



(x = 0,t) . 



Recursive application of the exponential operators with the aft-operator [H l , . . .]_ of the Cartan sub-algebra 
in (|3,14p specifies the form of the time- like Lax matrix 9Jt(x, t) by using the components a? of the root vectors 
a of corresponding ladder operators E± 



m{x, t) := H l ai(x, t) + E a b a (x, t) = 

= W (a^\x = 0,t) + fdy (d t xf\y,t))) + 



aGroot 



a 4 (^(x,t);%) = Oi(x,t) = Af ni \x = 0,t)+ / dy (d t xf\y,t)) ; 

Jo 

b a (Mx,t);kj) = b a (x,t) =Bg ni >(x = 0,t) x &p{ J* d£ \f\£,t) a*} 



(3.15) 

£ ex^j f d£ Af l) (e,0 «*} £ Q 4^(x = 0,t) ; 

Crnnt ^ 

(3.16) 
(3.17) 



We summarize the algebraic properties of the zero-curvature condition with following additional notes : 
If £(x, t) belongs to the generators of a closed algebra with the physical fields being the parameters or angles 
of the corresponding Lie group, one can generate 9Jt(x, t) for the zero-curvature condition by (|3.14l3.15p which 
then determines the nonlinear equations (|3.16|3.17p as GP-type equations, or other types corresponding to the 
chosen algebra of £(x, t) C sl(n, C). One has to take into account that one does not obtain an overdetermined 
system which finally results with the achieved matrix 9JT(x, t) f|3. 1413. 15[) and its chosen parametric dependence 
a,i(<j)p(x,t);kj), b a ((j)p(x,t);kj) f|3. 12|) into contradictory equations constraining the solutions from the zero- 
curvature condition to fixed time- and spatial-terms without a physical time development. 



3.2 Gauge invariance of the Lax pair and the zero-curvature condition 

It has already been stated that the Lax pair £(x,t), 9Jt(x,t) and its zero-curvature condition is by no means 
unique. If one restricts to transformations (|3. 18|) with a solely time dependent gauge matrix So = 9(t), one can 
immediately verify the invariance (|3.2U|) of the zero-curvature relation under the solely time dependent gauge 



3.2 Gauge invariance of the Lax pair and the zero-curvature condition 
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transformation (|3,19p 

So := S(i) ; {d x %) = {d x S(t) ) = ; (3.18) 
£^£' = 5 £% 1 ; 971 ^ 97t' = So S 1 + (^So) S 1 5 (3.19) 

(d t £) - (d x vri) + [£ , art] _ = o (d t £') - {d x m') + [£' , sW] _ = o . (3.20) 

The gauge invariance with So = 9(x,t) (|3.21II3.28~1) can also be extended to a general spacetime dependence 
within the corresponding Lie group, following from the closed algebra £(x, t), as e. g. in the Cartan-Weyl basis. 
As we begin with the transformations (|3. 23113 .25]) of So = S{x, t), acting onto the auxiliary field E(x, t) and Lax 
pair £(x,t), Wl(x,t), we attain the invariance of the Lax pair equations (|3.26p and its zero-curvature relation 
(13.271) . provided that the Maurer-Cartan relation (13.28P is fulfilled by the gauge matrix So = 9(x,t) 

Gauge invariance So = S(x, t) with Lie group of the closed algebra from £(x,t) ; 

(d x E) = £E; (d t E) =WIE; (3.21) 

=> (d t £) - (d x Wl) + [£ , Wl] _ = ; (3.22) 

E -» S' = S H; (3.23) 

£ - £' = So £ S 1 + {d x %) S 1 ; (3.24) 

art on' = So 9Jt s 1 + ($So) S 1 ; (3.25) 

=> (3*3') =£'»'; (d t H') =OJt'E'; (3.26) 

- (d x M) + [£' , SUf] _ = ; (3.27) 

provided that : c^So) Sq 1 ) - 0*(($So) Sq 1 ) + [(^So) So* , ($So) S = . (3.28) 

On the condition of the general Maurer-Cartan relation (I3.28|) . we can conclude for a whole set of equivalent 
Lax matrices (|3.23M3.2"6|) and their zero-curvature relations (|3.27p . However, it is in general even possible to 
transform the spatial Lax matrix £(x,t) with a suitable gauge matrix Q(x,t) to a diagonal form t) = 

H l Xf 1 \x,t) constrained to the commuting Cartan sub-algebra within a Cartan-Weyl basis of a Lie algebra. 
The diagonal form £^(x, t) of the Lax matrix £{x,t) has been used in section [37X1 to derive a general solution 
of the zero-curvature condition for the corresponding set of physical, nonlinear equations. In the following we 
construct the equations for the gauge matrix 9(x,t) and its algebra q 1 with 'rotation' angles i?j(x,i) which 
transform the general Lax matrix £(x,t) to a diagonal form £jj(x,t) = H l Xf 1 \x,t) 

£#M) = 9(x,t) £(x,t) T\x,t) + (d s S(x,t)) S~\x,t) ; (3.29) 
Q (x,t) = Q j ^(x,t); S(x,t) =exp{ (x,t)} . (3.30) 

In order to simplify the Lie group current (d x $(x,t) ) ^)~ l {x,t) in (|3.29p . we apply relation (|3.3ip which allows 
to resolve the current into the action of the aD-operator [Q(x,t) , ...]_ onto dg(x,t) / dx with the function 

(e x - l)/x 



{d x S(x,t)) S~W) = [d x exp{g(x,t)}j exp { -g(x,t)} = (3.31) 

f 1 i r t m dg(x,t) , . f 1 , r — i-n dg(x,t) 

= dv exp{vQ(x,t)} — — — exp { - v g(x,t)} = / dv exp [g(x, t) , . . . ]_ j — — — 

> 

( exp{[g(x,t) , ...]_}- 1 dQ(x,t) 

V [g(x,t), dx 
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After insertion of (|3.3ip into (|3.29p . we achieve eq. (|3.32p for the transformation of £(x,t) — > £ft(x,t) and 
dg(x,t) I dx. Under the assumption of an invertible function (|3.34|) for the aO-operator [q(x, t) , ...]-, one can 
finally derive a first order differential equation (j3.33j) for the Lie algebra g(x,t) = g 1 $j(x,t) of the gauge 
matrix $(x,t) so that the general spatial Lax matrix £(x,t) can be converted to a diagonal form £a(x, t) = 
JET* Xf 1 \x,t) with a suitable angular dependence $j(x,t) within the gauge transformation 



dg(x,t) 



(exp{[ (M) ,...]_} £(x,i)) + 



exp{[ (x,t) , ...]_}-! dg(x,t) 



dx 



dx 



1 



/exp{[fl(x,t) , ...]_}- 1 

V [g(x,t), ...]_ 

e x - 



-l 



£ 6 (x, t) - exp { [g(x, t) , ...]_} £(x, t) 



x 



(3.32) 

(3.33) 
(3.34) 



4 Determination and independence of conserved quantities 
4.1 Calculation of conserved quantities from the Lax pair 

The Lax pair £(x,i), SPt(x,t) and its zero-curvature condition, which specifies the nonlinear equations of 
physical fields, is accompanied by conserved, time-independent quantities. Let us consider the matrix T(x, t; k) 
J]) or the monodromy matrix T(i; k) (|4.2p with periodic, spatial boundary conditions on a circle x S [0, 2ir) 



%{x,t;k) = expj J d£ £(£, i; &)} ; 
1(t; k) = ^Ep{^ d£ £(£,£;&)} 



(4.1) 
(4.2) 



then the supposition, that all fields are periodic in x with period 2tt, implies that traces of powers of the 
monodromy matrix generate conserved quantities C^ n \k) = C^ n \t; k) independent of time 

C^(t;k) = Tr[(l(t;k)) n ] . (4.3) 

In order to attest this statement, we outline the defining, spatial ordering (j4.4H of the exponential with generator 
£(£, i; A;) of the spatial Lax matrix 

%(t;k) = expj J d£ £(£,i; &)} = exp{A*£(£ = 2n,t; k)} exp{Ax£(£ = 2vr - Ax,t;k)} . . . x (4.4) 
x ... exp {ai £(£ = 2 Ax, t; k)} exp {ax £(£ = Ax, t; k)} . 



The time-like derivative (<9jT(i; k) ) (|4.5p of the monodromy matrix (|4.2p involves the product rule with 
(dt£(y, t; k) ) according to the spatial ordering of the exponentials (|4.4p . After substitution of (dt£(y, t; k) ) by 
the zero-curvature condition, the integrand reduces to a total, spatial derivative d y (. . .) (|4.5p within the inte- 
gration boundaries y G [0, 2ir) of a circle. Hence, we can perform the spatial (fy-integration along the circle and 
acquire a commutator (|4.7p between the time-like Lax matrix 9Jt(x = 0, t; k) at the origin and the monodromy 
matrix T(i; A;), due to the presupposed periodicity (|4.6I) on a spatial circle 



(d t Z(t; k) 



^ dy exp{ J dfc £(&,*;*;)} (d t 2(y,t;k)) exp { jf #1 fc)} (4.5) 



4.2 Involution of conserved quantities and the classical x-matrix 
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exp{ / d£ 2 £(&,*;*;)} W{y,t;k) exp{ d& £(6, t; fc)} 



j/=2tt 
y=0 



= = 2vr, t; k) %{t; k) - X(t; fc) SDT(x = 0, i; fe) ; 

circle : =^ Tt(x = 2tt, i; k) = Wi{x = 0, t; k) ; 
(d t 1(t;k)) = [Wl(x = 0, t; fc) , T(t; fc)] _ . 



(4.6) 
(4.7) 



The conversion (|4.7f) of the time-like derivative (<9jT(i; k)) to a commutator allows to demonstrate the time 
independence of traces of arbitrary powers of the monodromy matrix %{t\ k) (|4.2p 



a 4 Tr 



(T(i; k) ) 



n Tr 



(3T(t;fc)) n l = 0; = 



, n— 1 



n Tr 



[£0l(a; = 0, i; fc) , X(t; k)] _ (%(t; k) 



,n-l 



C^Ht-k) = C^ n \k) . 



■0; 
(4.8) 



Therefore, the general time independence of C^ n \t; k) = Tr[(T(i; k) ) n ] is confirmed for arbitrary powers n G N 
so that the C^ n \k) = C^ n \t;k) (|4.8j) have to be regarded as the conserved quantities within a Liouville 
integrability of the corresponding nonlinear equations for the physical fields. Instead of the label 'n' for 
powers of %(t;k) (|4.2p . one can also conduct a power series expansion with the spectral parameters kj in 
order to generate conserved quantities from (orthogonal) polynomials of kj. Note that in the described cases 
of the (1+1) GP-equations for the attractive and repulsive interactions with the su(2) and sp(2,R) algebra, 
respectively, 



£su ( 2)(Q,?;fc) 



-ik —ip* 
t/j ik 



G 



£sp( 2 ,R)(Q,?;£0 



—k ip* 
tp k 



(k e 



(4.9) 



the spatial Lax matrices essentially stay unaltered as one introduces the external potential V(x,t). One can 
even choose a completely constant spectral parameter k(x,t) := ko so that the conserved quantities C^ n '(t; k) 
are not affected by the inclusion of an external potential V(x, t) into the GP-equations with a solely attractive 
('su(2)') or repulsive ('sp(2,R)') interaction. Therefore, the independence of conserved quantities C^ n \t; k) can 
be directly transferred from the well-known case of the integrable, (1+1) GP-equations without any external 
potential to the more general GP-types with an arbitrary external potential V(x,t). 



4.2 Involution of conserved quantities and the classical r-matrix 

In order to prove the independence of the derived, conserved quantities C^ n \t;k) of previous section [4TTI one 
has to verify the involution from the Poisson brackets of the canonical fields (I4,10p , These have to be taken 
into account because the conserved quantities C^ n \t;k) are generators for first order evolution equations with 
corresponding parameters t n (I4.1ip and have furthermore to fulfill the Jacobi identity (I4.12l4.13p . following 
from the Poisson brackets of canonical fields 4> a (x,t) = (Q a (x,t) , T a (x,t) ) 



= {C^(t;k),C^(t;k)} ; (m, n 2 G N) ; 



(4.10) 
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<} = {C^\t;k), ...}; (4.11) 



dt ni 

= {C^\t;k) ,{C^\t; k), ...}}+ (4.12) 
+ {C^(t;k),{... ,C^(t;k)}} + {... ,{C^(t;k),{C (n ^(t;k)}} ; 

<s=> = ( _g - — j - { {C( ni )(t; fc) , C ( " 2) (t; fe)} , . . . } . (4.13) 

\dt ni dt n2 dt n2 dt ni J s v / 



The involution of traces of powers of the monodromy matrix 

{C^ ni \t;k) ,C (n2 \t;k)} = 0, (4.14) 

is usually investigated by the so-called r-matrix approach where one assumes the validity of a 'fundamental 
Poisson bracket relation 7 (|4.15p between tensor products £l(xi, t; ki), £ 2 (x 2 , t; k 2 ) of spatial Lax matrices (cf. 
appendix A in [3] and chap. 2.5 in [TO] ) 



{£i(fi,t;A:i) f £ 2 (&,t;k 2 )} = [t 12 (h, k 2 ;t) , h) + £ 2 (x 2 ,t;k 2 )) _ 5(6 ; (4.15) 

112(^1,^2;*) = -t2i(fe2,fci;t) . 

The Poisson bracket {£1(6, t; &i) f £2(^2, ^ ^2)} (|4.15p is replaced by the commutator with the ti2(fci, k 2 ;t)- 
matrix where we further require the ' ultralocal' form with the spatial delta function 5(£i —£2)- (The 'ultralocal' 
condition is in general obtained from spatial Lax matrices which only depend on the physical fields Q a (x,t), 
7 a {x,t) without any derivatives of these. This condition is hence fulfilled for the (1+1) GP-equations.) As 
we consider more general types of X(x, y; t; k) matrices (|4.16p with initial coordinate y and end point x for the 
monodromy matrix (|4.2p instead of T(x; t; k) (|4.ip 

T(x,y;t;k) = ^ (4.16) 

we can start out from the general, tensorial Poisson bracket relation (|4.17|) whose right-hand side results from 
the validity of the Leibniz product rule for the Poisson bracket operation { . . . f . . . } 



{Xi(xi,yi;t;A:i) ? T 2 (x 2 , y 2 ; t; k 2 )} = / d& / d& Xi(xi,fr;t; A*) 1 2 (x 2 ,fr,t;k 2 ) x (4.17) 

Jyi Jy-2 

x {£1(6,*; **) ? £2(6,^^2)} Xi(Ciiyi;t;*i) ^2(6,1^;*;^) ■ 



In appendix[Bl we demonstrate according to Ref. [10] how to achieve the involution {C^ ni \t; k),C^ n2 ^ (t; k)} = 
from the generally valid relation (|4.17p of extended monodromy matrices %{x, y; t; k) under the special assump- 
tion of the 'fundamental Poisson bracket relation' (|4.15p with a spatially constant ti 2 (ki,k 2 ; i)-matrix. Since the 
spatial Lax matrices £ su (2)(Qi J 3 ! k), £ S p(2,R) (Q> k) of the integrable, (1+1) GP-equations do not change under 
the inclusion of an arbitrary external potential V(x, t) (apart from a possibly chosen spacetime dependence 
of the spectral parameter k — > k(x,t)), the whole derivation of the involution of conserved quantities can be 
directly conveyed from the case without an external potential to the case with an arbitrary potential V(x,t). 

The given derivation of the involution property (|4.10p of conserved quantities C^ (t; k) in appendix [B] 
depends on the assumed fundamental Poisson bracket relation (|4.15p with a spatially constant t\ 2 {ki,k 2 \t)- 
matrix. In the following we attain more general statements as we begin from relation (|4.17p . caused by the 
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Leibniz product rule of the tensorial Poisson brackets, and apply the gauge invariance of the Lax matrices 
which allows to transform to the diagonal, commuting Cartan sub-algebra. Further transformations of above 
relation (|4.17p rely on simplification of the tensorial Poisson bracket { . . . f . . . } of the spatial Lax matrices 
£i(x\, t; k\), £2(3:2, t; k 2 ). As we suppose a closed Lie algebra for £(x,t;k), we can diagonalize or parametrize 
to A(x, t; k) (of the commuting Cartan sub-algebra) with the invertible 'eigenvector' matrix HX(x, t; k) which 
also contains the ladder operators E±. This kind of parameters is taken within each part of the tensor product 
space £i(x\,t;ki), £ 2 (x 2 , t;k 2 ) of the spatial Lax matrix £(x, t;k) 

£(x,t;k) = ii(x,t; k) A(x,t;k) it" 1 (x, fc) ; A(x, t; k) = W Xf l) (x, t; k) ; (4.18) 
£i(a?i,t;A^) = Ui(a?i,t; h) A^x^t; fa) i^Vi) t; fa) ; (4.19) 
£ 2 {x 2 ,t;k 2 ) = ii 2 (x 2 ,t;k 2 ) A 2 (x 2 ,t;k 2 )ii 2 1 (x 2 ,t;k 2 ) . (4.20) 

On the condition of the "ultralocal" case of spatial Lax matrices, one can transfer the tensorial Poisson 
bracket (|4.21|) { . . . f . . . } of Lax matrices £1(2:1, t; fa), £ 2 (x 2 ,t;k 2 ) (I4.19)4.20p to terms with their eigen- 
values Ai(x%, t;k%), A 2 (x 2 ,t; k 2 ) and 'eigenvector' matrices Ui(xi, t; ki), il 2 (x 2 ,t; k 2 ) and a further part con- 
sisting of general spacetime dependent £12(2:1, fci; x 2 , k 2 ; t)-, t 2 \(x 2 , k 2 ; x\, fa; t)-matrices within commutators of 
£1(2:1, t; ki), £ 2 (x 2 ,t; k 2 ), respectively. The Xi 2 (xi,fa;x 2 ,k 2 ;t)-, t 2 i(x 2 , &2! 2:1, &i; i)-matrices f|4.22|) are com- 
posed of further parts &12, D12 (I4.23|4.24j) which contain terms with the Poisson brackets among the eigenvector 
matrices Hi(xi, t; fa), il 2 (x 2 ,t; k 2 ) and between the eigenvector matrix iii(x%, t; fa) and the diagonal Cartan 
sub-algebra part A 2 (x 2 ,t; k 2 ). In summary the basic nine terms in (|4.2ip are achieved by the Leibniz product 
rule of the tensorial Poisson bracket of canonical fields and can be grouped into a first part with the Pois- 
son bracket of the eigenvalues {Ai(xi, t; k\) f A 2 (x 2 ,t;k 2 )} and a second part with commutators between 
the X\ 2 {xi,fa;x 2 ,k 2 ;t)-, t2i(x2, k 2 ; x±, fa; t)-matrices and the corresponding spatial Lax matrices £i(x±,t; fa), 
£ 2 (x 2 ,t; k 2 ) (in the ultralocal case '5(xi — x 2 )') 

{£i( Xl ,t;fa) f £ 2 {x 2 ,t;k 2 )} = (4.21) 
= ii 1 (x 1 ,t;k 1 )iX 2 (x 2 ,t; k 2 ) {Ai(xi,t; fci) f A 2 (x 2 , t;k 2 )} ii^ 1 (x 1 , t;fa) ii 2 1 (x 2 , t;k 2 ) + 
+ ([ti2(x 1 ,ki;x 2 ,k 2 ;t) , £1(2:1, t; fa) ] _ - [x 2 i(x 2 , k 2 ; x\, fa; t) , £ 2 (x 2 , t; k 2 ) ] _) 5(xi - x 2 ) ; 

tx 2 {x u ki\x 2 M;t) = Oi2 + ^ [£12 , £2(2:2, t;k 2 ))_ ; (4.22) 
£12 = {Hi(xi,t;fci) f U 2 (x 2 ,t;k 2 )}i^ 1 {x 1 ,t;fa)iJ^ 1 {x 2 ,t;k 2 ) ; (4.23) 
O12 = iX 2 (x 2 ,t;fe2){ili(xi,t;A;i) f A 2 (x 2 ,t;k 2 )} ^{x^fa) ^ l {x 2 ,t;k 2 ) . (4.24) 

The (spatially constant !) r-matrix approach, given in appendix [Bl only allows to conclude for the involution 
of the conserved quantities C^ n \t;k), provided that the eigenvalues Ai(a;i, t; fa), A 2 (x 2 ,t; k 2 ) of the spatial 
Lax matrices are in involution (or their tensorial Poisson brackets vanish completely) (|4.25p . This can be 
accomplished for a symmetrical dependence on canonical fields <p a {x,t) = (Q a (x,t) , 7 a {x,t)) within the 

eigenvalues Xf 1 ) (Q a ,T a ; k) = X[ 6 \"? a ,Q a ;k) which may originate from symmetrically chosen parameters 
4> a (x,t) = (Q a (x,t) , y a (x,t) ) within the original Lax matrix £(Q a , Jq,; kj) of a closed algebra determined by 
a Cartan- Weyl basis of ladder operators (e. g. 4> a E a = Q a E" + T a E") 

{Ai(2;i, t; fa) f A 2 (x 2 ,t;k 2 )} = § ] => ^symmetric dependence on (4.25) 
canonical fields (p a = (Q a (x,t) , !P a (x,t) ) as e.g. X { " 3) (x,t) = J kj W + Q a ■ T a ) . 
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Under the assumption of a symmetrical dependence on the canonical fields Q a (x,t) , J' a (x,t) (|4.25p . the 
general tensorial Poisson bracket (|4.21|) reduces to the commutator (|4.26j) between the ti 2 (xi, &i;:c 2 , & 2 ; £)-, 
t 2 i(x 2 , fa; X\, fci;f)-matrices and the spatial Lax matrices £i(xi, t; fa), £2(^2)^5^2) 

{£i(xi,t;fci) ? £ 2 (x 2 ,t;£; 2 )} = (4.26) 

= ([ti 2 (xi,/si;x 2 ,/c 2 ;t) , £i(xi,t; - [r2i(iC2,fc2;5ci,fci;*) , £2(^2, *; fo)] _) <5(^i - x 2 ) . 

Relation (|4.26j) is similar to the already assumed fundamental Poisson bracket eq. fj4. 15|) : however, the given 
^12(^1, k±)X2, fa;t)-, t 2 i(x 2 , fa; xi, fa; £)-matrices have a priori no special symmetries, neither anti-symmetric 
nor symmetric. In order to be applicable for the derivation in appendix [Bj we transform by a commutator 
between the spatial Lax matrix £\{x\,t; fa), £ 2 (x2, i; fa) an d an additional symmetric matrix &12 = o" 2 l A4.27I - 
I4.29P which retains eq. (|4.26p invariant. This allows to assign the anti-symmetry to the transformed ti 2 = —X21 
matrix ()4.30|4.3ip so that the derivation of appendix [B] can be performed with the further restriction to a 
spatially constant Xi2(fa, fa; t) = — *2i (fa ? fa ; f) matrix ()4.32p 

Xvi{xi,fa;x2,fa;t) -> v 12 (xi,k 1 ;x2,k 2 ;t) + [a 12 (x 1 ,k 1 ;x 2 ,k 2 ;t) , 2 2 (x2,t;k 2 )]_ ; (4.27) 

x 2 i(x 2 ,k 2 ;x 1 ,k 1 ;t) -> r 2i (x 2 , k 2 ; x\, fa; t) + [a 2 i(x 2 , fa; x\, fa; t) , £i(xi, t; ki)] _ ; (4.28) 

o"i 2 (xi, fa; x 2 , fa; t) = 0" 2 i(iC 2 , fa; x\, fa; t) ; (symmetric matrix !) ; (4.29) 

ti 2 (til;#2;i) + [<7 12 (t|l;ti2;i) , £ 2 (x2,t;fa)]_ = — (r 2 i(tt2; JJ1; t) + [<7 2 i(tt2; , £ 1 (x 1 ,t;fa)]_y, (4.30) 
+ r 12 (ttl; (J2; t) + t 2 i(|J2; tfl; t) = - [a 12 ((ll; tf2; t) , £i(si, t; fcj) + £ 2 (x 2 , t; fa)] _ ; (4.31) 
Problem of spatial dependence of : ti 2 (a;i, fa;X2, fa;t) — > Xi2{fa,fa;t) . (4.32) 

Apart from the involution of conserved quantities according to the so-called ti 2 -matrix approach in appendix 
IB1 and [TO] , we also suggest a different proof which relies on a gauge transformation of the spatial Lax matrix 
£(x, t; k) to a completely diagonal form £^(x, t; k) = A(x, t; k) given in section [331 and (j4.38|4.39|) . This results 
in vanishing n 2 -, t 2 i-matrices ()4.33|4.34p within the basic tensorial Poisson bracket relation (|4.21M.24"1) so 
that the Poisson bracket of the diagonal Cartan sub-algebra matrices £#,1(2;, t; k) = K\(x, t; k), £# )2 (x, t; k) = 
A 2 (x,i;/c) only remain because the diagonalizing matrices ili(x 2 , t; fa), lt 2 (xi, t; fa) reduce to tensorial unit 
matrices (|4.35p . We remark again that £(x,t;k) and its gauge transformed, diagonal form £f,(x,t;k) = 
A(x,t;k) represent the same nonlinear equations of fields 4> a (x,t) = (Q a (x,t) , T a (x,t)) within the general 
sl(ra, C) algebra or within one of its sub-algebras as e. g. su(n) 

if £{x,t;k) — > £fi(x, t; k) = A(x, t; k) by a gauge transformation (4.33) 

=> ti 2 (tJl;tJ2;t) =0 A r 2 i(tt2; |J1; t) = ; (4.34) 

{£i(x 1: t;fa) f £ 2 (x 2 ,t;A:2)} => (4.35) 
={Ai(x 1 ,t;fc 1 ) f A 2 (x 2 ,t;k2)} 

, * V 

=> {£ 6>1 ( Xl ,t;fa) f £ A , 2 (x 2 ,t;& 2 )} + 

+ ([ti 2 (|jl;|j2;t) , £ 6tl (x u t; fa) ] _ - [r 21 (tJ2; J}l;t) , £ A , 2 (x 2 , t; fa) ] _) <5(xi - x 2 ) . 
=0 =0 
As we repeat to choose symmetrical dependences (|4,36p on canonical fields <p a (x,t) = (Q a (x,t) , T a (x,t) ) with 
diagonal, traceless Cartan sub-algebra matrices £jj[x,t; k) = H % \\ H ) (Q a ,'J > ol ; k) = A(Q Q , "? a ; k) 

\\ H \Q a (x,t) , y a (x,t)) =^> symmetrical dependences on Q a (x,t) , T a (x,t) ! ; (4.36) 
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\f\Q a (x,t), 7 a (x,t)) 



X\" >(y a (x,t) , Q a (x,t)); 



(4.37) 



we resolve the involution of the conserved quantities following from (|4.21l4.17p by the remaining, vanishing 
tensorial Poisson bracket (|4.37|) of the diagonal Lax matrices £^ i(Q a , 3^; k), £jj 2 (Q Q , y Q ; A:). In section I3T21 it 
has already been exemplified that the gauge transformation (|4.38l4,39p , which constrains the spatial Lax matrix 
£(x,t;k) to the Cartan sub-algebra £ft(x,t;k) = H % \\ H \Q a ,'J > a ;k), should exist for very general conditions 
of the physical fields 



dg(x,t) 
dx 



5(x,t) 



/exp{[g(x,t) , 
^ [g(x,t), 



-l 



H l A- , {Q a (x,t),7 a (x,t)) + 



exp{[g(x,t) ,...]_} (W ki + E«? a (x,t) + E° Q a (x,t) 
exp{g(x,t)} gauge transformstion should exist ! . 



(4.38) 



(4.39) 



5 Extension of the zero-curvature condition beyond (l+l)-dimensions 

5.1 Determination of the nonlinear equations for the fields of relevant, physical observ- 
ables 

In previous sections [2][5] we have emphasized the algebraic properties of Lax pairs in (1+1) dimensions and have 
also considered cases where one can obtain chaotic behaviour within the general, non-compact sl(n, C) algebra 
by separating into sub-algebra parts as e. g. su(n) or sp(n,M), etc. . In this section we point out a possible 
extension to ((N-l)+l) Euclidean dimensions x u = (x,t), x u = (x,t) for Lax matrices W(x,t) taking values 
within the sl(n, C)- or one of its sub-algebras. Instead of two equations for the (1+1) dimensions, one has to 
regard ./V equations (15, ip with Lax matrices W(x,t) acting onto the n-component, auxiliary field H(x, t). In 
analogy to a N dimensional electromagnetic theory, the Lax matrices W{x,t) are termed as matrix-potentials 
W(x, t) which depend on canonical fields <t> u a {x, t) = (Q^(x, t) , "J"^(x, t) ) and spectral parameters k\ (x, t) within 
a Cartan- Weyl basis (|5.2|3.8p of the general sl(n, C) algebra or a chosen sub-algebra as e. g. su(ra) 

8 u E(x,t) = W(x,t)E(x,t); (k, \, fi, v = 0, (1, . . . , AT — 1) ) ; (5.1) 
= ^ i ^(x,t) + J E Q <(x,t) = J ff J ^(x,t)+^^(f,t)+^Q^(f,t) ; (5.2) 
Matrix-potentials 21" = 21" (Q, IP; k) of canonical fields <j) v a {x, t) = {Q? a {x, t; k) , 7 v a (x y t;k)) . 

In comparison to previous sections [2][2J we assign to the time- like matrix potential 21 (x, t) the Lax matrix 
Wt(x, t) and to the (N-l) spatial matrix potential components 2l l (x, t), '2l(x, t)' the spatial Lax "vector" £(x,t) 
with corresponding (N-l) spatial Lax matrix components £*(x, t) (|5.3p . The equivalence of mixed and exchanged 

partial derivatives d tJ- d u E(x,t) = d u d fl E(x,t) of (|5.ip enforces N(N — l)/2 zero-curvature equations (|5.4p 
instead of the single zero-curvature condition in previous sections for (1+1) dimensions. In analogy to the 
electrodynamic case, one can introduce a matrix-valued field strength tensor ^ v and a matrix-valued current 
(15. 5p . given by the commutator of two matrix potentials 2P, 21" 



2l°(:z,i) = Wi{x,t) , W{x,t) = £ i (x,t) ; (* = 1, . . . , N - 1) ; 

= &*W - d v W + [21" , 2P] _ ; SI 1 ' G sub-algebra C sl(n, C) ; 



(5.3) 
(5.4) 
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5 EXTENSION OF THE ZERO-CURVATURE CONDITION BEYOND (1+1)-DIMENSI0NS 



= ffiyu - d v W = T u = , 21"] _ . (5.5) 

By using the Maurer-Cartan structure equations for the gauge matrices So 

(0"(d"S o ) So 1 ) " (0"(0"So) So" 1 ) + l[d v %) So 1 , (^So) So 1 ]. = , (5.6) 

we can straightforwardly derive a gauge invariance of the zero-curvature eqs. (|5.4)5.5p under the gauge trans- 
formation (|5.7p for the matrix potentials W — * %' v 

W - 21'" = So % v So 1 + (a"S ) S 1 ; (5.7) 
= - d v %'^ + [21'" , 21'^] _ ; So G sub-group C SL(n, C) . (5.8) 

In later sections we have to perform a gauge transformation in order to prove the involution of conserved 
quantities. One has to take a gauge transformation for the scalar product £(x(£),t) = (dx(Q I <K) • ^(x(C),t) 
along chosen 'loops' x(( = 0) = x(C = 2tt) to the Cartan sub-algebra elements £#(:?(£),£) = H l Xf 1 \x((),t) ■ 
(dx(()/dO. 

The nonlinear equations from (15.3ll5.5p can be transferred to the structure constants i c? , of the underlying 
algebra from the matrix potentials 2P(x, t) which are defined with the generators q j and corresponding fields 
^(Q, T; k) (]5.9p as an ansatz for nonlinear equations. As we resolve relations (|5.4|5.10p in terms of the generator 

basis $ and structure constants i <^ k l (j5.9f) . we finally attain the set of nonlinear equations (|5,lip with the 
ansatz of the ^(Q, k) fields which can also be related to the more fundamental fields Q^(x,t), J"^(x,t) in 
eq. ([521) 

[ d j ,S> k ]_ = i^q 1 ; W(x,t)=g i *$(Q t 9 ] k); (5.9) 
= (sf((d»*fi - {d^f)) + [tf , g k ]_ *J *jA ; (5.10) 
= (3 M ^(Q, T; k)) - (d v $?(Q, T; k)) + % c ik l #J(Q, IP; k) *£(Q, 5*; k) . (5.11) 



5.2 Conserved quantities of monodromy matrix paths with nontrivial homotopy of fields 

In the following we consider a closed 'loop' x(£) = x%, £ G [0, 2ir) (|5.12p within the group manifold of the chosen 
Lie algebra (sl(n, C) or a corresponding sub-group) for the matrix potentials 21" = (SEJl, £) where a reference 
point xp defines the beginning and the end of the loop 

x(0 = x e ; £ G [0, 2vr) ; x(£ = 0) = x Po ; x(£ = 2vr) = x Po . (5.12) 

It is of crucial importance that the fibre space, given by the N-l spatial coordinates for the base space and 
by the mapping to the group manifold of the generators 21* = £*, does not allow a continuous contraction of 
the loop to a trivial, single point as e. g. the reference point xp . This supposition yields a straightforward 
derivation of conserved quantities along nontrivial, non-contractable loops in analogy to section FPl One begins 
with the matrix path from the reference point xp to a point x(£) = x^, which is specified by a parametrization 
x((), (C G [0, £), £ G [0, 2-7r)), and performs a spatial ordering of exponentials exp{AC (dx^ / d() ■ £(x^,t;k)} 
along a part of the loop 



T(x£,x Po ;t;k) = 




(5.13) 



5.3 Involution of conserved quantities and the classical x-matrix in arbitrary spacetime dimensions 
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In correspondence to section 14. H we take the time derivative of f)5. 13|) and have to regard the product rule 
following from the spatial ordering of the exponential step operators along the part ( € [0, £) of the loop 

x c = f(C) 



d t %(x£,x Po ;t;k) = J d(exp^J 



dCi 



(5.14) 



~d( 



(d t £(x t;k))) / d£ 



C 



We can replace the scalar product (dx^ / d() ■ (dt£(x£,t; k) ) in (|5.14j) by the zero-curvature relations (|5.15|5.16|) 



dt£ - dWl = [an , £ 

dx^ 



d( 



•3*£ 



(5.15) 
(5.16) 



in order to transform the integrand of (|5.14p to a total derivative of the loop parameter £ € [0, £) (cf. section 
Oeqs. (I4.5ll4.7p ) 



d t 1(x^,xp )t;k) 



d£ exp 



r/.r 



Ci <v^ 



£{x Cl ,t;k) \ x 



(5.17) 



da?£ 
"df 



an,-i • £ 



cxp 



dC2 



dx 



<2 A 



d(2 



dC |(elq5 



x 9Jt(x£,t;/c) exp 



C 



d(i ^.£(%,t;fc) 



C 



dC 2 



dCi 
d% 
dC 2 



£(%,£; fc) 



= 3Jt(z(£),i; fc) 1(x^,x Po ]t;k) - %(x^,x Po ;t;k) Wl(x Po ,t;k) . 

As one takes £ = 2-7T so that x(£ = 2tt) = x(£ = 0) = xp are the beginning and the end of a closed loop '0' in 
the base space for a nontrivial, (non-contractable) homotopic mapping of the fields, one transforms the time 
derivative of monodromy matrix (|5.19p to a commutator which finally gives the conserved quantities (|5.2U|) of 
traces of powers of the monodromy matrix T(Q, t; k) 



%(x{£ = 2n),xp ;t;k) 
d t 1(Q,t;k) 

d t Tr[Z n (Q,t;k)] 



1(Q,t;k) ; 

[Tl(x Po ,t;k) ,%(Q,t;k)]_; 

{d t C^(Q,t;k)) =nTr\(d t 1{®,t;k)) 1 n_1 (0,t; k) 



(5.18) 
(5.19) 

(5.20) 



n Tr 



[m(x Po ,t;k) , 1(0, t; k)] _ 1(0, t; k) 



5.3 Involution of conserved quantities and the classical r-matrix in arbitrary spacetime 
dimensions 

As we have generalized the derivation of conserved quantities from a circle x G [0, 27r) within (1+1) dimensions 
to a closed loop x(£), £ € [0, 27r) within a (N-l) dimensional base space for a nontrivial, homotopic mapping of 
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6 SUMMARY AND CONCLUSION 



a Lie group SL(n, C) (or a sub-group as SU(n), etc. ), we can also extend the involution properties of section [4T2l 
for the (1+1) dimensional case to those of the ((N-l)+l) spacetime. One can repeat the calculations of section 
14.21 by replacing the spatial Lax matrix £(x,t) by the scalar product H(x^,t;k) = ■ £(x^,t;k) where the 
'loop' parameter £ € [0, 2ir) substitutes the spatial coordinate 'x' or '£' of the circle for the (1+1) dimensional 
case in section 14^21 

. — (%X/- -* 

Apply £(x£, t; k) = ■ £(x^,t; k) instead of £(x,t;k) or £(£, t;k) in section 14.21 and (5.21) 

repeat calculation of involution properties ; 
Diagonalize by a gauge transformation : (5.22) 

£(x c ,t;k) - 2 A (x t;k) = H i \f\x c ,t;k)-^ . 

Similarly, the classical r-matrix approach can be conveyed from appendix [B] of the (1+1) dimensional case to 
((N-l)+l) dimensions with the replacement (|5.2ip under restriction of spatially constant matrices x(k^,kj]t). 
However, we emphasize again that the analogous, corresponding transformations in place of the (1+1) dimen- 
sional case of section l4~2l and appendix [B] only hold for non-contr actable, nontrivial, homotopic mappings from 
the loop within the (N-l) dimensional base space to the chosen Lie group manifold as SL(n, C) or one of its 
sub-groups. 



6 Summary and conclusion 

6.1 Lax pairs and chaotic behaviour of (1+1) GP-type equations 

This article has been initiated by the notion whether any Lax pair construction can only lead to a completely 
integrable behaviour. As we have verified in sections 12.11 and 12.21 for Lax pairs of the (1+1) GP-equations as 
generators of the sl(2, C) algebra, one can even determine Lax pairs for arbitrary external potentials V(x, t) 
without changing the spatial Lax matrix component £(x,i). Since the conserved quantities and their invo- 
lution only depend on the exponential step operators with the spatial Lax matrix £(£,t), £ E [0, 2ir), one 
can directly conclude for a Liouville integrability as in the cases without an external potential, either from a 
spatially constant ti2-matrix approach according to appendix [B] or from a gauge transformation to the diago- 
nal, commuting Cartan sub-algebra elements as the eigenvalues of £(£, t) within the tensorial Poisson bracket 
relation (A symmetric dependence of the eigenvalues on the canonical fields within the Poisson bracket has 
to be presupposed.). As we reduce the generators of £(x,t), Wl(x,i), either to the sub-algebra su(2) or to 
sl(2,R) (of the most general, nontrivial sl(2,C) algebra, cf. appendix [A]) for an attractive or repulsive in- 
teraction, the hermitian property of the prevailing Hamiltonian, following from su(2) or sl(2,R) Lax pairs, 
prevents any chaotic behaviour, due to the chosen compactness with a spatial circle x S [0, 2ir). However, 
as we combine the ip S u(2)( x , t) an d ip S p(2,«.)( x > t) fields of the two integrable, (1+1) GP-equations with attrac- 
tive and repulsive interaction to the complex-valued parameter fields within the most general, non-compact 
sl(2,C) algebra, probability or density of the ip gu (2)(x, *) an d ipsp(2,R){ x i t) fields can flow and change between 
the coupled GP-equations which separately contain incoherent, non-hermitian terms for a chaotic behaviour. 
This chaotic behaviour from Lax pair construction is even possible for sl(n > 2, C) algebras where one has 
to select a sub-algebra for a compact sub-group as SU(n) C SL(n, C) so that coupled nonlinear equations of 
physical fields are also composed of incoherent terms, giving rise to unlimited increase of the ip su ( n )(x,t) fields 
and corresponding chaotic behaviour. 



6.2 Lax pair construction in arbitrary spacetime 



25 



6.2 Lax pair construction in arbitrary spacetime 

The given construction of Lax pairs for the (1+1) GP-equations as generators of sl(n, C) (or of a sub-algebra 
as su(n)) straightforwardly generalize to arbitrary spacetime dimensions. However, we emphasize again that 
this extension beyond (1+1) dimensions necessarily has to involve a nontrivial homotopic mapping from the 
loop within the (N-l) dimensional base space to the considered group manifold, following from the Lax pair 
generators. In absence of a nontrivial homotopy, it is possible to contract the loop within the fibre space 
to trivial point mappings so that the construction of the conserved quantities and their involution becomes 
trivial and meaningless. Therefore, Lax pair constructions for a Liouville integrability or a possible chaotic 
behaviour beyond (1+1) dimensions have to be accompanied by an investigation for a nontrivial homotopy of 
the underlying fibre space (|16j). 



A Reduction of gl(n, C) to sl(n, C) Lax pairs by separating the trivial 
trace parts 

In this part [A] of the appendix we assume that Lax matrices £(x,t), 9Jt(x,t) are not traceless and therefore 
belong to the gl(n,C) algebra as the most general case of n x n matrices. The general Lax matrices £(x,t), 
Wl(x,t) are separated into diagonal unity parts 1 A£(x, t), 1 AWl(x, t) and remaining traceless parts £o(x,£), 
DJto(x, t) (|A,l)A,2h of the sl(n, C) algebra (n = N% = N^yi). This trace separation is also performed for the initial 
matrix 9Jtj n j(x = 0,i) at the coordinate origin within the solution of the zero-curvature relation for 9Jt(x,t) = 
Wlo(x,t) + 1 ATl(x,t) where explicit use is made for the trace splitting of £(x,i) = £o(x,t) + 1 A£(x,t) with 
the 'oD '-operator 



£(x,t) = 2, (x,t) + j^-Tr[£(x,t)]=£o(x,t) + iA£(x,t); [Tr[£ (x,t)] = Oj ; (A.l) 

m(x,t) = m Q {x,t) + ^ ■Tr[Tl(x,t)] =m Q (x,t) + 1 AdJl(x,t) ; (Tr[mo(x,t)] = o) ; (A.2) 

m ini (x = o,t) = mf n \(x = o, t) + i Am ini (x = o,t) ; (rr[iB$%(x = o, t)) = o) ; (A.3) 

m(x,t) = fm (x,t) + i A<m(x,t) = (aa) 

= e^{ J d£ [£ (£,i) + l A£(£,t), ...]-} (^(x = 0, t) + 1 A9Jt mi (x = 0, tfj + 

+ J^dyt^[J d£[£ (£,i) + lA£(£,t), ...]_} ((d t £ (y,t)) + 1 {d t A£(y,t) 



As we regard the complete vanishing of the 'a0 '-operator part [1 A£(£,i) , . . .]_ in (|A.4j) . one can conclude for 
the separation of the solution of the zero-curvature relation into the two independent parts (1A.5|A.6|) where 
equation (|A.5j) only consists of the total traceless sl(n, C) generators £o(x,i), 9Jto(a;, t), 5DT^(x = 0,t), and 
where equation (|A.6|) separately has the remaining fields A£(x, t), A9Jl(x, t), A9Jti n i(x = 0, t) from the diagonal 
unity part '1' 

M (x,t) = e3q5{ £ d£ [£ (£, t) , . . . ]-}m ( fJ i (x = 0, t) + (A.5) 
dy^p[j dt[£o(Z,t), ...]_} [d t £ (y,t)) ; 

AWl(x,t) = Am ini (x = 0,t) + / dy (d t A£(y,t)) (A.6) 

Jo 
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B INVOLUTION OF MONODROMY MATRICES FOR SPATIALLY CONSTANT Vf\-MATRIX 



AMini(x = 0, t) + & ( dy Tr [£(y, t)] /N £ 



We can also verify from the zero-curvature relation (|A.7p the separation property into traceless sl(n, C) matrices 
and diagonal unit parts. Since the diagonal unit parts A£(x,t), A9Jt(x,t) only cause vanishing commutators 
within the general zero-curvature condition, the latter relation splits into the solely traceless part of sl(n, C) 
generators £o(x,t), 9Ro(x,t) and a simple part for the completely diagonal terms with the fields A£(x,t), 

Am{x,t) 

=0 



(dt£>) - {d x Wl) + [£ , 971] _ = (d t & ) - (d x dJl ) + [£o , St«o] _ + 

i 



(A.7) 



+ 



Nz 



8 t (Tr [£(x, *)] )-d x d t (j* dy Tr [Z(y, t 







(N £ = N m ) 



Therefore, the purely traceless generators with £o(x,t) (and the derived or traceless chosen generators for 
dJlo(x,t)) can give rise to nontrivial, nonlinear equations within a relevant Lax pair construction. 



B Involution of monodromy matrices for spatially constant t-matrix 

A standard proof of conserved quantities from the spatial Lax matrix £(x,t) C sl(n,C) begins with the 
assumption of the fundamental, tensorial Poisson bracket relation (jB.ip . having a delta function 5{x\ — X2) for 
the 'ultralocal' case of fields. We note that this relation transforms the quadratic term of £(x, t) to a linear part 
within a commutator of the so-called \\ 2 '-matrix which is supposed to be independent on the space coordinate 

V 

f 2 2 {x2,t;k 2 )} = [n 2 (fci, k 2 ;t) , £i(xi,i; h) + £ 2 {x 2 ,t; k 2 )] _ S( Xl - x 2 ) . (B.l) 

We consider the spatial evolution eqs. (|B,3)B,4j) for the initial coordinate 'y' and end point 'x' of the defined 
type (|B.2p of monodromy matrix T(x, y; t; k) which follows from subsequent spatial ordering of exponential 
step operators with Lax matrix £(£, t; k). As one applies the Leibniz and product rule for the tensorial Poisson 
bracket (|B.5j) of two monodromy matrices (|B.2[) . we accomplish the fundamental Poisson bracket (1B.1|B.6|) 
within a double spatial integral £1 € [yi,xi], £2 S [y 2 ,x 2 ] 

%{x,y;t-k) = &p{J de£(e,t;fe)} ; (B.2) 

(d x 1(x,y;t;k)) = £(x,t; k) T(x,y;t; k) ; (B.3) 
{d y 1(x, y;t;k)) = - %{x, y; t; k) £(y, t; k) ; (B.4) 

{Xi(xi,yi;t; h) ? T 2 (x 2 , y 2 ; t; k 2 )} = / d£i / d£ 2 %\{xu£y;t',kx) T 2 (x 2 , £ 2 ; t; fc 2 ) x (B.5) 



x {fiifo,*;**) ? £ 2 (6,i;fc 2 )} £2(6,1*2;*; fc 2 ) ; 

{£i(6,*;fci) f £2(6,^2)} = [ti2(*;i,fc2;t),£i(Ci,t;fci) + £2(6,t;fc2)]_*(ei-C2). (B.6) 

After substitution and insertion of (|B.1|B,6|B,3|B,4|) into (|B.5|) . the tensorial Poisson bracket of monodromy 
matrices reduces to simple integrals ([B.7IB.8P because the corresponding integrands only consist of total deriva- 
tives <5(£i — £2) (<% + %,) and d%, respectively 

{%l(x 1 ,y 1 ;t;k 1 ) f T 2 (x 2 , y 2 ; t; k 2 )} = %i(x 1 , x min ; t; h) T 2 (x 2 , x min ; t; k 2 ) x (B.7) 
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d£i ^2 5(6 - £ 2 ) 



2-1 {Xmim £ 2 ( 



x ti 2 (^i, & 2 ;i) f(% ^1(^1, y maa; ;t;fci)) + (%T 2 (6,ymaa!;i;fc2))J + 

+ ((%Ti(x mi „,^i;t; fei) ) + (%1 2 (smin,6;i; fa) )) ti 2 (fci, fe;*) x 

x £i(£i, Umax; t; h) 12(^2, Umax] t;k 2 ) 1i(y m ax,yi;t;ki) I2{y m ax,y2;t;k2) 

= t Xi(x 1 ,x min ;t;k 1 ) l2(x 2 ,x m in;t;k 2 ) / d£% d£, 2 S(Ci - 6) (% + %) 



£i]t;ki) T 2 (» m i n ,C 2 ;*; A; 2 ) ti 2 (fci, /c 2 ; Ti(^, y maa; ; t; fci) T 2 (£ 2 , y maa; ; fc 2 ) 
x 1i(y m ax,yi;t; fa) 1 2 (y max 1 

y2]t;k 2 ) ; 

{^(a;!,^;*; fci) f 2: 2 (x 2 ,y 2 ;i;fc 2 )} = (B.8) 
= Xi(xi,x min ;i; fci) c l 2 (x 2 ,x min ;t;k 2 ) / a?£ x 

x 9 € ^Ti(x min ,^;t;/ci) T 2 (x min , £; i; fc 2 ) ti 2 (fci, & 2 ; i) Tifoymaa;; fci) 1 2 (£, y ma x] t; k 2 

X 'XliVmax, yi) t'i k\) T 2 (y 

mux ! 

After performing the spatial integration in (|B.8p . we achieve the tensorial Poisson bracket (|B.9[) of evoluton 
matrices with the remaining integration boundaries for the spatial £ integration. Let us assume the integer 
relation (m > n 2 ) so that taking the integration boundaries transforms the right-hand side of (|B.9p to a 
commutator-like relation (|B,10p with the classical 'ti 2 '-matrix 



{Ti(2vrni,0;t;fci) f T 2 (2tt n 2 , 0; t; k 2 )} = 

= Xi(27rni,^;t;fci) 1 2 {2tt n 2 ,£;t;k 2 ) x 12 {k 1 ,k 2 ;t) £i(£,0;t;fci) T 2 (£, 0; i; /c 2 ; 

(5: 1 (27rm,0;t;fci) ? T 2 (27r n 2 , 0; t; k 2 )} = 
= r £i(2-Kni,2-Kn2;t;ki) x 12 (ki, k 2 ; t) Ti(27r n 2 , 0; t; k\) T 2 (27r n 2 , 0; t; k 2 ) + 
- 1i(27rni,0;t; fci) X 2 (2-7r n 2 , 0; i; k 2 ) t\ 2 {ki, k 2 ;t) . 



£= 2 7r Min(n\ ,712) 



(B.9) 



(B.10) 



As we perform the total trace Tri 2 in order to obtain the conserved quantities C^ ni \k\), C^ n2 ^(/c 2 ) within the 
Poisson bracket of physical, canonical fields, the right-hand side of (IB.lOj) yields the trace TVi 2 of a commutator 
with the 'ti2 '-matrix which results into zero and therefore demonstrates the independence of the conserved 
quantities C( ni )(A;i), C {n2 \k 2 ) 



Tr 



12 



{T 1 (27rn 1 ,0;t;A: 1 ) ? T 2 (2vr n 2 , 0; t; k 2 )}\ = {C^(h) , C^\k 2 )} = 

Ii(27rni,0;*;fci) 

\' * ' 

Tri 2 Ti(27rn 2 ,0;t; fci) %\{2-k m,2ir n 2 ;t; k\) 1 2 (2TTn 2 ,0;t;k 2 ) ti 2 (k\, k 2 ; t) + 

Ti(27rni,0;t; k\) X 2 (27r n 2 , 0; t; k 2 ) ti 2 (fci, k 2 ; t) = 



(B.ll) 



Tr 
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Xi(27rrti,0;t; k±) T 2 (27r n 2 , 0; t; k 2 ) , ti 2 (&i, & 2 ; i) 



= . 



28 



REFERENCES 



The given proof of this appendix iBl relies on the spatial independence of the 'ti2 '-matrix, ti2(&i, fe; t) with 
the ultralocal property 5{x\ — x^) in the assumed fundamental Poisson bracket relation (|B.llB.6p , The spatial 
independence of \\i '-matrix has to be incorporated in order to transform the spatial integrations in (|B.7IB.8P 
to total spatial derivatives so that the integrands simplify to remaining integration boundaries in (|B.9IB.10p . 
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